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Basic  models  for  electromagnetic  transients  arising  in  exoatmospheric  pho¬ 
ton  bursts  were  developed  25  years  ago  [1,2].  Explicit  numerical  calcula¬ 
tions  have  been  carried  out  since  then  [3,4].  These  calculations,  based  on  a 
"nominal"  (i.e.,  monoenergetic)  photon  source,  predict  rather  large  field 
intensities,  even  at  the  ground  level.  We  find,  however,  that  this  result  is 
not  to  be  expected  if  we  assume  a  spectrally  finite  photon  source.  In  typical 
applications,  a  spectrally  finite  case  is  more  likely  [5,  p  36].  One  might  pre¬ 
sume  that,  even  in  such  a  case,  the  monochromatic  source  is  the  proper 
one,  since  kinetic  processes  involving  the  electromagnetic  field  sources 
(x-ray  photoproduction,  Compton  scattering,  electron  positron  pair  pro¬ 
duction)  are  quantum  processes  and  may  therefore  be  taken  to  be  incoher¬ 
ent.  Nevertheless,  the  time  scale  for  these  processes  (<10-12  s)  is  much 
smaller  than  the  transient  pulse  time  scale  (~10-8  s),  so  the  fields  may  be 
assumed  to  originate  from  coherent  sources. 

The  problem  under  consideration  may  be  envisioned  if  we  start  with  a 
wave  equation  for  the  transverse  fields  in  the  impulse  approximation  [6]: 

9[rE(r,r)]_-Z0rj(r,T) 

dr  2  '  '  ' 

where  t  is  the  retarded  time  [f  -  (r/ c)]  and  Z0  is  /JqC,  the  impedance  of  free 
space.  The  photon  pulse  gives  rise  to  a  primary  electron  current,  which  de¬ 
pends  on  the  kinetic  process  for  electron  production.  In  turn,  the  primary 
current  generates  a  secondary  current  through  the  process  of  ionization  by 
collision  with  ambient  molecules.  The  magnitude  of  the  secondary  current, 
which  is  generally  proportional  to  the  primary  current,  will  normally  build 
up  to  the  magnitude  of  the  primary  current,  so  that  the  field  intensity  from 
that  point  onward  propagates  as  a  free  wave  whose  magnitude  is  essen¬ 
tially  independent  of  the  source.  This  point,  then,  is  referred  to  as  the  satu¬ 
ration  limit.  We  find  this  condition  schematically  by  setting  the  ratio  of  the 
magnitude  of  the  primary  current  to  that  of  the  secondary  current  equal  to 
unity.  If  we  use  the  expression  ;  =  -enpvp  for  the  primary  current,  and  the 
expression^  =  -enf>s  for  the  secondary  current,  we  have  for  the  saturation 
condition: 


(2) 


The  first  term  on  the  left  is  obtained  from  the  well-known  formula  for  ion¬ 
ization  by  collision  [7]: 


ife) 

d%  [dnp  j 


=m , 


(3) 


where  £  is  a  path  length,  and  T(j8)  is  a  function  of  =  (vp/ c)  and  the  other 
parameters  given  by  Fermi  [7].  We  then  have 
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nJL=  [to  Y{p)pcdT 

nP  Jo  l-P 


(4) 


We  integrate  over  the  time  profile  of  the  incident  pulse  (in  retarded  time), 
and  use  for  r  a  value  equal  to  tQ  =  10-8  s.  Thus,  equation  (2)  becomes 


r  wQ 

l-p 


vs(E)  =  1  , 


(5) 


and  the  order  of  magnitude  of  the  saturated  field  is  determined  as  soon  as 
we  establish  the  relation  between  the  drift  velocity  and  the  field  intensity. 

Note  that  this  result  is  essentially  independent  of  the  intensity  of  the 
source  and  that  the  role  of  the  geomagnetic  field  is  of  secondary  impor¬ 
tance.  Only  in  the  case  of  x-ray  photoelectrons  do  we  have  to  consider  the 
details  of  the  synchrotron  motion  in  the  impulse  limit.  For  relativistic  elec¬ 
trons,  the  motion  in  the  impulse  limit  is  nearly  a  straight  line,  and  the  role 
of  the  geomagnetic  field  is  less  significant  in  the  sense  that  it  is  the  deter¬ 
mining  factor  for  the  directional  effects,  but  not  the  field  intensity.  Since 
the  most  intense  fields  are  produced  by  the  relativistic  electrons,  we  will 
ignore  the  effect  of  synchrotron  motion  in  our  discussion. 


The  relation  between  drift  velocity  and  field  intensity  needed  in  equa¬ 
tion  (5)  is  rather  complicated.  The  simplest  possibility  that  can  be  used  as  a 
basis  for  discussion  is  a  linear  mobility  model: 

vs  =  ^E  >  (6) 

so  that,  finally,  from  equation  (5)  we  obtain  for  the  saturation  field 


p .  1  zl 

nmtQ  •  (7) 

The  value  for  mobility  of  6  m2/V-s  is  suggested  by  Longmire  et  al  [4].  It  is 
found  to  be  surprisingly  good  for  the  conditions  under  consideration  be¬ 
cause  mobility  is  expected  to  be  both  field-  and  time-dependent.  The  drift 
velocities  as  normally  reported  are  obtained  using  quiescent  conditions, 
and  thus  may  not  be  suitable  for  use  in  transient  calculations  [8].  The  pro¬ 
cedure  used  in  the  numerical  work  reported  by  Longmire  et  al  [4]  was  that 
of  calculating  the  mobility  at  each  stage,  and  updating  its  value  in  each  it¬ 
eration.  Detailed  numerical  calculations  so  obtained  agreed  in  order  of 
magnitude  with  the  numerical  estimates  found  by  using  a  fixed  value  for 
mobility.  We  will  use  the  value  of  6  m2/V-s,  as  suggested,  in  making  our 
own  estimates,  and  justify  this  approximation  from  another  viewpoint 
later  on. 


Our  contention  is  that  the  procedure  of  using  a  monoenergetic  photon 
source  leads  to  an  overestimate  of  the  field  intensities  to  be  expected  at  low 
altitudes.  Photons  that  propagate  down  toward  the  earth  encounter  an 
increasingly  dense  array  of  target  molecules  with  molecular  density  n(z)  at 


r 
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altitude  z,  determined  by  the  Law  of  Atmospheres  for  the  regions  of 
interest  [9], 

n(z)  =  n( 0)  e<  ,  (8) 

where  £  =  (z/A),  with  A  a  parameter  whose  value  is  about  7  km,  and  n( 0)  is 
the  molecular  density  at  ground  level.  As  photons  propagate  into  the 
lower  atmosphere,  reaction  rates  increase  and  the  photon  flux  decreases. 
Ultimately  a  peak  reaction  level  is  reached,  where  reaction  rates  are 
maximum  [10]: 


e 


(9) 


where  M.w)  is  the  mean  free  path  for  collisions  at  ground  level,  i.e.,  A(w)  = 
[l/(n(0)o(ro))].  (We  omit  directionally  dependent  effects  in  the  above  rela¬ 
tion  by  restricting  ourselves  to  observations  of  pulses  that  occur  directly 
overhead.)  The  energy-dependent  cross  section  is  equal  to  the  Thompson 
cross  section,  cT,  multiplied  by  a  parameter  G(w)  with  w  =  (j hv/mc 2),  i.e., 
o(w)  =  OjG{w)  [11],  so  that 

e^  =  G(w)e^T  .  (10) 

Here,  is  obtained  from  equation  (8)  using  the  Thompson  cross  section. 
The  numerical  value  of  obtained  this  way  is  5.2.  This  corresponds  to  a 
height  of  around  36  km.  G(w)  is  a  monotonically  decreasing  function  of  w. 
For  x-ray  photoproduction  G(w)  >  1,  and  the  reaction  zone  ranges  from  an 
altitude  of  around  80  km  to  about  36  km.  For  Compton  scattering  G(w)  <  1, 
and  the  reaction  zone  ranges  from  about  36  km  to  about  18  km.  The  high- 
intensity  fields  created  at  altitudes  of  30  to  36  km  are  attenuated  by  the  sec¬ 
ondary  currents  created  by  the  high-energy  electrons  at  altitudes  of  18  to 
24  km.  This  can  be  seen  by  referring  to  equation  (7)  and  noting  the  /)  =  1  for 
ultrarelativistic  electrons  that  are  produced  at  these  altitudes. 

We  proceeded  in  our  estimate  of  the  field  intensities  by  replacing  the 
monoenergetic  source  with  a  distributed  one  characterized  by  a  spectral 
density  n(w)  such  that  jdw  n(w )  =  1,  with  n(w)  determined  from  the  data  in 
figure  1  [5,  p  361].  Current  density  in  the  reaction  zone  was  calculated  by 
making  a  saddle-point  expansion  of  variables  around  the  height  defined  in 
equation  (9)  and  using  spectral  averages,  i.e.,;'  =  jdw  n(w)jp(w)  and js  = 
jdw  n(w)js(w),  etc.  We  obtained  primary  current  densities  that  were  fairly 
uniform  over  the  range  of  altitudes  in  the  reaction  zone,  but  the  secondary 
currents  in  the  reaction  zone  decrease  rapidly  with  altitude.  This  leads  to 
an  attenuation  of  the  field  intensity  as  it  propagates  down  to  the  Earth's 
surface. 


Numerically  computed  saturation  field  levels  directly  below  the  source  lo¬ 
cation  (linear  mobility  model)  are  presented  in  figure  2.  The  calculations 
were  made  using  a  mobility  of  6  m2/ V-s,  and  the  results  for  the  ground  in¬ 
tensities  were  found  to  be  severely  attenuated  by  the  secondary  electrons 
associated  with  the  high-energy  tail  of  the  photon  spectrum,  i.e.,  for  pho¬ 
tons  above  8  MeV.  Our  numerical  estimate  for  the  ground-level  field  inten¬ 
sity  is  about  6300  V/m. 
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Figure  1.  Spectrum  of 
initial  photon  energy 
2000  yd  from  20-kT 
explosion  [5,  p  361]. 
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Figure  2.  Saturated 
electric  field  versus 
photon  energy  for 
linear  and  nonlinear 
models.  Field  value 
for  particular  photon 
energy  value 
represents  magni¬ 
tude  of  saturated 
electric  field  that 
would  be  found  if 
photon  spectral 
distribution  were 
truncated  for  all 
larger  photon  energy 
values. 
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Figure  3.  Comparison 
of  saturated  electric 
field  predicted  by 
our  linear  model  and 
Longmire  model  for 
1.6-MeV  mono- 
energetic  gamma 
source.  Longmire 
data  obtained  from 
Longmire  [3,  fig.  13]. 


Apart  from  including  the  effects  of  the  complete  photon  spectrum,  our 
model  does  not  differ  significantly  from  the  Longmire  model  [3].  By  using 
a  limiting  procedure,  our  method  recovers  the  monochromatic  results  for 
an  impulse  source  of  1.6  MeV  gamma  energy.  A  comparison  of  the  pre¬ 
dicted  saturation  field  level  is  presented  in  figure  3  for  both  models  (data 
obtained  from  Longmire  [3],  fig.  13). 

An  improved  estimate  of  the  swarm  parameters  is  possible  if  one  consid¬ 
ers  the  pseudo  transport  model. 


dv  _  eE  v 
dt  m  tc  ' 


(ID 


where  tc  is  a  collision  time.  If,  for  tc  we  take  (X/v),  we  find  for  the  steady- 
state  drift  velocity  the  well-known  result  [12], 


v  =  CE  (1+a)/2  , 


(12) 


where  a  represents  some  velocity  dependence  of  the  collision  cross  sec¬ 
tion.  Since  C  is  density-dependent,  one  finds 


(l+a)/2 


(13) 


with  K  a  constant  and  p  in  atmospheres,  to  be  a  better  representation.  Ulti¬ 
mately  one  has,  empirically,  for  oxygen,  a  good  fit  for  the  drift  velocity, 
given  by 


(14) 
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From  this  we  can  calculate  field  mobilities  and  we  find  the  value  of 
6  m2/ V-s  to  be  a  good  approximation  for  the  range  of  values  in  the  region 
of  interest.  We  can  also  use  these  parameters  to  obtain  an  estimate  of  the 
time  for  equilibrium  to  set  in  for  the  plasma.  We  find  this  to  be  less  than 
10-9  s,  except  for  electrons  of  lowest  energy.  It  is  thus  possible  to  go  back  to 
equation  (1)  and  use  equation  (14)  to  calculate  the  saturation  field  (with,  of 
course,  spectrally  averaged  quantities).  We  can  set 

hec  =  ~ensec^  '  (15) 

with  X  =  1400V  G(w)  in  MKS  units. 

Since  the  wave  equation  is  now  nonlinear  and  not  even  separable,  this  is  a 
major  task.  However,  along  geomagnetic  flux  lines,  the  equation  is  ap¬ 
proximately  separable.  From  equation  (2),  we  develop  a  model  wave  equa¬ 
tion  of  the  form 


e  H(1  -  E)  =  exp  -  J  .  (18) 

We  calculate  ra  from  equation  (17)  using  spectral  averages  and  we  find  it 
to  be  about  7.8  km.  Typically  (r;  -  r)  =  20  km,  which  makes  the  right-hand 
side  e~ 2,6  or  =  0.08,  and  E  =  1  in  equation  (16).  By  the  use  of  spectral  aver¬ 
ages  we  find  in  equation  (16)  that  £j  =  7000  V/m. 

One  sees  here  that  the  saturation  condition  E  =  1  is  not  fulfilled.  As  a  mat¬ 
ter  of  fact,  E  ranges  in  value  from  about  0.5  to  1  over  the  range  of  param¬ 
eters  of  interest.  Nevertheless,  figure  2  indicates  that  this  model  provides 
consistent  results  with  the  linear  mobility  model  when  the  source  includes 
a  high-energy  tail. 

We  have  modified  the  standard  model  [3,4]  for  predicting  the  electric  field 
intensity  from  an  exoatmospheric  gamma  source  by  including  the  effects 
of  all  components  of  the  photon  pulse.  We  have  found  that  the  high- 
energy  components  of  the  photon  pulse,  which  constitute  only  a  small 
fraction  of  the  total  spectrum,  are  responsible  for  a  reduction  in  the  field  at 
ground  level. 
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Attn  4200  Ctr  for  Sp  Sensing 

Attn  4600  Condensed  Matter  &  Radiation  Sci 

Attn  4700  Plasma  Physics  Div 

Attn  5700  Tactical  Elect  Warfare  Div 

Attn  6800  Elect  Sci  &  Techlgy  Div 

Attn  8000  Ctr  for  Sp  Techlgy 

Attn  Code  4820  Techl  Info  Div 

4555  Overlook  Ave  SW 

Washington  DC  20375-5000 

Cmdr 

Nav  Surfc  Weapons  Ctr 
Attn  Code  E23 1  Techl  Lib 
Dahlgren  VA  22448-5020 

Hqdtrs  USAFA/DFSELD 
Attn  HQUSAFA/DFSELD 
2354  Fairchild  Dr  Ste  3A22 
USAF  Academy  CO  80840-6214 

US  Air  Force  Phillips  Lab 
Attn  PL/WSR  Baum 
3550  Aberdeen  Ave  SE 
Albuquerque  NM  87116-008 
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Distribution 


Lawrence  Livermore  Natl  Lab 
Attn  L-156  M  Bland 
AttnL-86  HS  Cabayan 
PO  Box  808 
Livermore  CA  94550 

Sandia  Natl  Lab 

Attn  Orgn  3141  Reports  Acqstn 

PO  Box  5800 

Albuquerque  NM  87185 

Natl  Inst  of  Stand  &  Techlgy 
Attn  V  Ulbrecht  Rsrch  Info  Ctr 
RmEOl  Bldg  101 
Gaithersburg  MD  20899 

West  Virginia  Univ 
Dept  of  Physics 
Attn  DR  AD  Levine 
Morgantown  WV  26506 

Fed  Commctn  Comm 
Attn  DR  A  H  Rausch 
2000  M  Stret 
Washington  DC  20554 

US  Army  Rsrch  Laboratory 
Attn  AMSRL-SL-C  S  Share 
Aberdeen  MD  21005 

US  Army  Rsrch  Lab 

Attn  AMSRL-SL-C  W  Hughes 

Aberdeen  Proving  Ground  MD  21005-5066 

US  Army  Rsrch  Lab 

Attn  AMSRL-SL-CO  A  Van  de  Wal 

Aberdeen  Proving  Ground  MD  21005-5067 


US  Army  Rsrch  Lab 

Attn  AMSRL-SL-CM  M  Mar 

Aberdeen  Proving  Ground  MD  21005-5068 

US  Army  Rsrch  Lab 

Attn  AMSRL-WT-TB  R  Lottero 

Aberdeen  Proving  Ground  MD  21005-5425 

US  Army  Rsrch  Lab 
Attn  AMSRL-SL-CS  J  Beilfuss 
Attn  AMSRL-SL-CM  D  Farenwald 
Edgewood  Arsenal  MD  21010-5423 

US  Army  Rsrch  Lab 

Attn  AMSRL-OP-SD-TA  Mail  &  Records 
Mgmt 

Attn  AMSRL-OP-SD-TL  Tech  Library 
(3  copies) 

Attn  AMSRL-OP-SD-TP  Tech  Pub  (5  copies) 

Attn  AMSRL-PS-PD  J  McGarrity 

Attn  AMSRL-PS-PD  T  Oldham 

Attn  AMSRL-WT-ND  B  Luu 

Attn  AMSRL-WT-ND  C  S  Kenyon 

Attn  AMSRL-WT-ND  G  Merkel 

Attn  AMSRL-WT-ND  J  R  Miletta 

Attn  AMSRL-WT-ND  R  J  Chase  (20  copies) 

Attn  AMSRL-WT-NF  L  Jasper 

Attn  AMSRL-WT-NH  R  A  Kehs 

Adelphi  MD  20783-1 197 
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